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Abstract — In this work, we consider the problem of detecting the 
presence of a new user in a direct-sequence/code-division-multiple- 
access (DS/CDMA) system with a doubly-dispersive fading channel, 
and we propose a novel blind detection strategy which only requires 
knowledge of the spreading code of the user to be detected, but no 
prior information as to the time-varying channel impulse response and 
the structure of the multiaccess interference. The proposed detector 
has a bounded constant false alarm rate (CFAR) under the design 
assumptions, while providing satisfactory detection performance even 
in the presence of strong cochannel interference and high user mobility. 

Index Terms — DS/CDMA, doubly-dispersive channel, user detection, 
cell-search, neighbor discovery. 

I. Introduction 

Code-division multiple-access (CDMA) is widely employed in 
third-generation (3G) cellular standards, and is a strong candidate 
for the definition of the air-interface of future wireless networks 
[l]-[3]. In this paper, we investigate the problem of detecting the 
presence of a new active user which naturally arises in the design 
of CDMA-based systems. Indeed, it is known that hidden (i.e., 
undetected) active terminals can cause a substantial degradation 
of the bit-error-rate performance even if the optimum maximum- 
likelihood multiuser receiver is employed [4]. If users enter and/or 
leave the channel at random epochs due to the bursty nature of 
the data transmission, detecting and exploiting these changes in the 
user population is crucial to achieve high spectral efficiencies and 
to guarantee a minimum quality of service to the set of active users. 
The problem of user acquisition may also be reframed in the context 
of cell-search and handoff procedures [5], wherein mobiles have to 
detect strong surrounding base stations and have to synchronize with 
their unique signatures transmitted on the broadcast channel, and in 
the context of neighbor discovery [6], [7], wherein mobiles sense 
the surrounding environment to identify other nodes. 

First works on the detection of a spread-spectrum signal embed- 
ded in a mixture of multiaccess interference and thermal noise have 
focused on synchronous CDMA systems operating over non-fading 
or frequency-flat fading channels. In [8] a two-stage procedure is 
proposed, wherein the first stage resorts to linear or non-linear fil- 
tering to remove the multiple-access interference, while the second 
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stage performs detection. In [9] a generalized likelihood detector 
is proposed, while a detection procedure based on the generalized 
cumulative sum test is presented in [10]. Finally, [11] proposes 
an innovative data detection structure based on the random set 
theory: despite its good performance, this approach has a prohibitive 
computational complexity. The work in [12], instead, discuss strate- 
gies for user detection in asynchronous DS/CDMA systems with a 
frequency-selective, slow fading channel: the proposed procedures 
are based on the application of the generalized likelihood ratio test 
(GLRT). 

All these previous works critically rely on the assumption that 
the channel conditions remain constant for many symbol intervals. 
However, this assumption may not be met in practice since user 
mobility and terrain changes may induce rapid variations in the 
channel response, which become more pronounced as the carrier 
frequency increases [13] — [15]. The problem of achieving reliable 
communication over rapidly time-varying channels naturally arises 
for example in many popular data transmission standards (e.g., 
W-CDMA, UMTS, DVB-T and IEEE802.16), and in vehicular 
networks that have become more and more popular lately. Hence, 
it is of interest to develop improved algorithms that are capable of 
operating in this more hostile environment. Many papers tackling the 
problem of reliable communication over doubly-dispersive channels 
have appeared in the literature, [ 16]— [ 1 8] just to cite the most recent 
ones. However, the majority of these previous studies focus on 
orthogonal frequency division multiplexing (OFDM) systems, and 
do not address the problem of new user detection. 

In this work, we consider a DS/CDMA system operating over a 
doubly-dispersive fading channel, and we assume that neither pilot 
symbol nor secondary data are available to estimate the channel 
impulse response and the structure of the cochannel interference. 
Under this challenging scenario, we first show that the GLRT 
approach followed in [12] cannot be employed anymore. Instead, 
we resort to the so-called method of sieves, a valuable statistical 
tool reported in [19], [20], and we derive a novel modified GLRT 
(MGLRT) which is fully-blind (i.e., it only requires knowledge of 
the spreading code of the user to be detected) and has a bounded 
constant false alarm rate (CFAR) under the design assumptions (i.e., 
the test statistic, although not enjoying itself the CFAR property, is 
upper-bounded by a CFAR test). Simulation results indicate that 
the proposed MGLRT is robust to the presence of strong cochannel 
interference and performs well even when the user to be detected 
is not active during the entire processing interval. 

The remainder of the paper is organized as follows. In Section II 
the system model is presented. Section III contains the synthesis of 
the new detector. Section IV contains the numerical results, while 
concluding remarks are given in Section V. 

Notation: In the following, X T , X\ X + and rank{X} denote 
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the transpose, the Hermitian transpose, the Moore-Penrose pseu- 
doinverse, and the rank of a matrix X G C mxn , respectively. 
tr{A}, | A | and \A\ P denote the trace, the determinant, and the 
product of the positive eigenvalues of a square matrix A, respec- 
tively. I m and On.m indicate the identity matrix of order m, and 
an n x m-dimensional matrix with all-zero entries, respectively. 

II. System model 

Consider a DS/CDMA ^system with K asynchronous users and 
denote by s k (t) = J2 n =o Pk(n)iptx(t - nT c ) the spreading 
signature assigned to user fc, wherein N is the processing gain, 
T c is the chip interval, k = [f3 k {0), . . .,(3 k (N - 1)] T G is 
the spreading code, and ^tx(i) is a unit-energy pulse waveform 
that is non-zero in [0, PT C ) with P a positive integer. The complex 
envelope of the received signal can be written as 



A'-l +oo 



!/(*) 



k — q= — oo 
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x c k (t, r)dr + w(t) . 



Tfc 
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In (1), r k > is the transmission delay of the fc— th user; 
b k (q) 6 B k is the data symbol transmitted by the fc— th user 
during the q— th symbol interval of length Tt = NT C , where B k 
is an arbitrary (possibly user dependent) finite constellation with 
average unit energy; c k (t, r) is the equivalent time- variant channel 
response of the fc-th user (accounting for the transmit signal power, 
the path loss and the multipath fading); finally, w(t) is the additive 
noise, modeled as a circularly-symmetric, zero-mean white Gaussian 
process with power spectral density Mo- 

The received signal y(t) is first sent to a linear filter matched to 
iptx{t). Letting T m be the maximum multipath delay spread of the 



channel, ip(t) = 



V'tx(w)V'tx( tt — t + PT c )du, and assuming 



at the design stage that c k (t,r) remains constant for about Tt, + 
T m + PT C seconds, the output of the receive filter can be written 

as 

r(t) = / y(u)ip^(u - t + PT c )du 



-1 +oo 



(2) 



= X1 E b k (q)J2^(n)g k (t-nT c ;q)+n(t). 



fc=0 q= 



wherein n(t) is the filtered noise process, and we have shoved the 
channel-dependent quantities of the fc-th user into the function 

/+oo 
ij)(t-T-Tk- qT b )c k (qT b + r k , r)dr . (3) 
-oo 

Notice that g k (t; q) has compact support in [qTb + r k , qTb + r k + 
T m + 2PT C [ C [qT b , (q + l)T b + 2PT C ), where the inclusion on 
the right hand side stems from the fairly reasonable assumption 1 

T k + Tm < T b . 

After sampling r(t) at rate T c /M, with M a positive integer 
representing the number of samples per chip, we stack the complex- 
valued samples corresponding to the interval [qTb, (q + L)Tb), with 
L > 2, into the LNM— dimensional data vector r(q). Skipping the 

'This assumption can be easily relaxed as shown in [21], [22]. 



analytical details, r(q) can be expressed as 

K-l L-l 

r(q) = EE ! " (9 + £)Ck,l9k{q + £) + n(q) , (4) 

k=0 l=-2 

wherein g k (q) G (£(n+2P)m . g a c i ian ije[.j e p e i 1 dent vector which 
contains the samples g k (qTb + nT c /M; q) for n = 0, . . . , (N + 
2P)M — 1; n(q) G <C LNM is a complex circularly-symmetric Gaus- 
sian vector with zero-mean and full-rank covariance matrix R n G 
C iJVMxiJVM , with {R n ) i3 = M i>((i-j)T c /M + PT c ); 
{C M G C LJVMx(JV+2P)M ', t = L - 1} are code- 

dependant matrices which contain suitable shifts of the spreading 



vector (3 k . Indeed, let A k 
matrix 



/3 fe (0) 
: /3 fc (0) 
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be the 
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' Im , (5) 
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the matrix Ck,o is obtained by taking the first LNM rows of A k , 
while the matrices {C k: e, I ^ 0} are obtained by taking suitable 
sub-blocks of Ck.o [12], [21]. It is worth pointing out that the 
discrete-time signal representation (4) is very powerful since it 
isolates the known code-dependant matrices from the unknown data 
symbols and the unknown channel-dependent vectors. 

III. Detector design 

In order to cope with the time-varying nature of the channel, 
Q consecutive data vectors are jointly processed by the detector. 
Thus, the problem of detecting the presence of a new user, say user 
0, in the received data multiplex can be formulated in terms of the 
following binary hypothesis test: 



H :r(q) = w{q), q = 1, 

Hi : r(q) = b (q)C , a g a (q) + z{q), q = 1, 



,Q, 



(6) 



where 



w(q) = ^2 ^ b k (q + £)C k , e g k (q + I) + n{q) , 

fc=l l=-2 

L-l 

z(q) = w(q)+ ^ bo(q + £)Ca,iga(q + I) . 

t=-2 

Henceforth we assume Q > LNM, which implies that the data 
matrix R = [r(l), . . . , r(Q)] G C LNMxQ is full row rank 
with probability one under both hypotheses. Also, under Ho, we 
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model w(l), . . . ,w(Q) as independent and identically distributed 
(i.i.d.) complex circularly-symmetric Gaussian vectors with zero- 
mean and covariance matrix M w , which is unknown since no 
prior information as to the interferers number and their channel 
impulse response is assumed available; similarly, under Hi, we 
model 2(1), . . . , z(Q) as i.i.d. complex circularly-symmetric Gaus- 
sian vectors with zero-mean and unknown covariance matrix M z . 
Finally, we assume that only the spreading code /3o of the user to 
be detected is known, while its data symbols {foo(l), ■ ■ ■ ,bo(Q)} 
and its channel-dependent vectors {go(l), . . . , go(Q)} are regarded 
as unknown deterministic parameters. 

Lacking prior information as to the user parameters and the 
second-order statistics of the interference, the Neyman-Pearson test 
[23] cannot be implemented here. A common way to circumvent 
this drawback is resorting to a GLRT approach which amounts 
to replacing the unknown quantities with their maximum likeli- 
hood (ML) estimates under each hypothesis [23]. Upon defining 
W = [w{l)...w(Q)] 6 C LMNxQ , Z = [z(l)...z(Q)] G 
C iMArx Q and q = [bo( i) 9o (i) 5 . . . , bo(Q)go(Q)} G C DxQ with 
D = (N + 2P)M, the test (6) can be recast as 



H ■■ R = W 
Hi:R= Co.oG + Z 



and the GLRT can be expressed as 



sup f R (R\M x ,G,H!) Hl 



sup f R (R\M w , Ho) 



H a 



(7) 



(8) 



where r\ is a threshold to be set so as to achieve a given probability 
of false alarm Pf a and 

f R {R\M w ,H a ) = 



f-'o 



- 7r e X p{-tr{M+RR t }}d(^lR) , 



(9) 



m (R\M w ,H ) 



f R (R\M z ,G,Hi) = 

— exp {-tr {M+ (R - CG) (R - CG) f } } 



m 1 (R\M- e ,G,H 1 ) 



S (®l(R-CG)) 



are the probability density functions (pdf's) of the received data 
under both hypotheses [24]. In the above expressions, Co and ci are 
positive normalization constants; S(-) is the product of the Dirac 
delta functions of the elements of the matrix argument; finally, & x 
is a LNM x (LNM — rank{ M x }) matrix whose columns form 
an orthonormal basis for the null space of the columns of M x , with 
x — w, z. 

Unfortunately, the GLRT (8) cannot be applied in this scenario; in 
fact, after maximizing the numerator with respect to M z , we have: 



sup f R (R\M Zl G,Hi) 
oc sup 



(ID 



gT \(R-CG)(R-CG)% 



Now, the matrix G can always be chosen such that (R — CG)(R — 
CGy has D strictly positive eigenvalues approaching zero, thus 
causing the right hand side of (11) to diverge. As a result, the 
GLRT (8) is unbounded, and hence the ML estimate does not exist. 
This circumstance was also observed in [20, Equation (4)], [25, 
Equation (5)], and is explained by noticing that the parameter space 
is too large. A possible solution to this drawback is to resort to the 
method of sieves [19], wherein the parameter space is restricted to a 
subspace such that the ML estimate exists and is unique. According 
to the method of sieves, a unique ML solution under Hi exists upon 
restricting the parameter space to the set 



5 = 




= rank(M z ) < LNM - D 

B Op t 2NM-p 
02NM-p,p 02NM-p,2NM-p 



where M z = (R - CG)(R - CG)\ B is an arbitrary p x p 
matrix and \I' is the matrix of the normalized eigenvectors of M z , 
and we consider the following modified GLRT (MGLRT): 



sup mi(R\M z ,G,Hi) H 

(M„,G)eS > 

sup m (R\M w ,H ) < V ' 



(12) 



The next proposition provides a closed-form expression for the test 
statistic in (12). 

Proposition: The MGLRT for the binary detection problem in (6) 
is given by 



T = 



\RBl\ 



Hi 

\{Ilnm - CC+)RRi(I LNM - CC+)t|„ < ' 

-Ho 



(13) 



Proof: Using [20][Lemma 1], it follows that the test (14) can be 
equivalently recast as 



\R&\ >\ 

G : p<LNM-D \(R ~ CG)(R — CG)^\p < 

no 

Hence, it is sufficient to show that the solution to the problem 
G = arg inf \(R - CG)(R - CG) f L 

G : p<LNM-D 

is G — C + R. Let C = U'SV' 1 be the economy-size singular 
value decomposition (SVD) of C, where U G C LNMxD is a 
tall matrix with orthonormal columns, S G C DxD is a diagonal 
matrix containing the non-zero singular values on its main diagonal, 
and V G C DxD is a unitary matrix. Let F = . . . , f D ] H = 
Sy+G. Since SV f is invertible, we have that G = VE^F with 

F = arg min \<R - UF)(R - UF^L . 

F : p<LN RI — D 

Generalizing [25] [Proposition 1], it is verified that F — [/* R. This 
concludes the proof. 



Some remarks are now in order. Notice first that the projector 
I lnm — CC + spans the null space of C and, hence, the matrix 
{Ilnm — CC + )R in the denominator of (13) does not contain 
the useful signal CG. Also, the computation of the test statistic 
T can be simplified as outlined in the following. Let R = U'R, 



2 For notational simplicity, in the following we drop the subscript in Co,o- where U G 



iLNMxLNM 



is a unitary transformation such that 
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[/At/ 1 " and 
A 



Olnm-d,lnm-d Olnai-d,d 
Od,lnm-d Id 



Also, let [L, Olnm,q~lnm]Q be the Q-R decomposition of R 
with L G c LNMxLNM a lower triangular matrix and Q G C QxQ 
a unitary matrix. We have: 

\RR?\ 



r ■ 



' \(Ilnm- A)RRi(I LNM 
\(Ilnm - A)LL1(I LNM - 

LNM 

n (m 2 - 

i=LNM-D+l 



(14) 



where denotes the (i, i)-th entry of L. Notice now that com- 
puting the Q-R decomposition of R only involves 0(Q(LNM) 2 ) 
floating point operations [26]. 

A. On the detector's CFARness 

In the following, we investigate the ability of the proposed 
detector to operate with a CFAR under the design assumptions. 
In particular, let M* G £ LNMxLNM be the true covariance 
matrix of the received data under Ho and let PP^ be the Cholesky 
decomposition of UM*U f , with P G C iJVMxiJVA/ a lower trian- 
gular matrix. Then, under Ho, we can write l-RP^I = \P\ 2 \ VV'| 
and {I LNM - A)R = [V?P},Oq,d]\ where V G C LNMxQ 
is a random matrix whose elements are i.i.d. complex circularly- 
symmetric Gaussian variables with zero-mean and unit variance, 
Vi G C ( - LNM ~ D)xQ contains the first LNM - D rows of V and 
Pi G c (LNM - D)x( - LNM - D) contains the first LNM - D rows 
and columns of P. Finally, letting [L, Olnm,q-lnm]Q be the 
Q-R decomposition of V , the test statistic in (14) can be recast as 

\P\ 2 \VV^\ 



T 



\Pi\ 2 \ViV?\ 

2NM 

n 

i=LNM-D + l 



2 
Pi.i 



n 



(15) 



j = LNM-D+l 



where pij and h,i denote the (i,i)-th entry of P and L, respec- 
tively. 

Since T e depends upon P (and, hence, upon ), the proposed 
receiver does not have a CFAR. However, it is possible to make 
the receiver bounded-CFAR with respect to M^, through a suit- 
able normalization of the test statistic. Indeed, let T etTnax be the 
maximum value of T e over all possible scenarios under Ho (which 
can be obtained via experimental measurements). The normalized 
test statistic T n — T /T etmax is upper bounded under Ho by the 
CFAR test statistic 7cfar = T /T e . Therefore, if we choose r\ so 
as P(7cfar > T}\Ho) = Pfa, then we have 



P(T„ > T]\H ) < P(T C FAR > V\H0) = Pfa ■ 



IV. Numerical Results 



(16) 



The performance of the proposed detection scheme can be 
expressed in terms of the probability of false alarm Pf a and the 
probability of detection Pj. Since no closed form expressions are 



available for both Pf a and P^, we resort here to Monte-Carlo 
simulations. In the following, we consider a DS/CDMA system with 
oversampling factor M = 2, processing gain N — 15, pseudo-noise 
(PN) spreading codes and bandlimited raised cosine chip waveforms 
with roll-off factor a and duration 4T C . Each user transmits binary 
data symbols, i.e., Bu = { — 1, +1} V k, while the processed data 
vectors r(l), . . . , r(Q) span L = 2 symbol intervals with Q — 120. 
The equivalent baseband channel responses are modeled as 3-path 
channels with time- varying path gains [13], [14]: 



c k (t, t) = A k a ktP (t)8(T - T k ,p), 



,K. (17) 



In (17), A k denotes the received signal amplitude of the fc-th user; 
{•7~fc. p } denote the path delays and are modeled as independent 
random variables uniformly distributed in [0, (N — 1)T C ]; finally, 
{a k ,p(t)} denote the path gains and are modeled as independent 
stationary Gaussian processes with zero mean and autocorrelation 
function set according to the Jakes model [27], i.e. R(t) = 
F>[a ktP (t)a* k p (t - t)] = J (2-k fdT^r) , where J (-) is the zero 
order Bessel Function of the first kind and fd is the normalized 
Doppler frequency shift. In the following, we set fd to 0.1 or 
0.01, which correspond to high user mobility scenarios. We define 
the Signal-to-Noise-Ratio (SNR) as QA^/No', also, we assume 
A\ = ... = Ak and define the Signal-to-Interference-Ratio (SIR) 
as A 2 / A 2 . Unless otherwise stated, we consider a power-controlled 
scenario with SIR=0 dB and we set the test threshold so as to have 

Pfa = ID" 2 . 

In Fig. 1 we plot P& versus SNR for different values of K and 
a = 0.3. As a benchmark, we also report the detection performance 
of the genie-GLRT detector which has perfect knowledge of the 
covariance matrices M z and M w , i.e., 



sap f R (R\M X , G, H! 

G > 



f R (R\M w ,H ) 



(18) 



Ho 



Taking the logarithm of (18) and after some algebra, we obtain 

Q 

J2 r H<l) (M- 1 - M' 1 ) r(q) 

9=1 

Q Hi 
+ YV(g)M" 1 C (C^M^Cy 1 &M^r{q) > r? 

9 = 1 Ho 



(19) 



It is interesting to notice that the detection performance improves 
for larger fd since temporal diversity is exploited. Also, the per- 
formance loss of the MGLRT with respect to the genie-GLRT is 
less than 3 dB for K = 1 at any SNR, while it becomes larger for 
K > 1 as a consequence of the fact that the estimation accuracy 
of the unknown covariance matrices under Ho and Hi inevitably 
degrades. In Fig. 2, instead, we plot Pd versus SNR for different 
values of the roll-off factor and K — 3. As expected, the detection 
performance improves for increasing a since a larger bandwidth is 
employed. 

So far we have assumed that the user to be detected is active dur- 
ing the entire processing window. However, in real world situations, 
a new user may start transmitting asynchronously with respect to 
the beginning of the test. In particular, let Q < Q and assume that 
only the data vectors r(Q — Q + 1), . . . , r(Q) contain the useful 
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signal to be detected. Fig. 3 shows that the detection performance 
gracefully degrades for decreasing values of Q. This is a very 
attractive feature, which makes the MGLRT detector suitable for 
the detection of bursty users. 

Finally, in Fig. 4 we remove the power-controlled assumption 
and we study the detection performance of the MGLRT detector 
for different values of the SIR. It is seen that Pd is not significantly 
affected by the power level of the interfering signals. In fact, for 
negative SIRs the detector suffers no asymptotical loss with respect 
to the power-controlled scenario, whereas Pd improves for positive 
SIRs. As a consequence, the MGLRT detector can be effectively 
employed in near-far scenarios, as for example in neighbor discov- 
ery applications wherein no power control can usually be guaranteed 
among mobiles [7]. 

V. Conclusions 

In this paper we have considered the problem of detecting the 
presence of a new user in a DS/CDMA system operating over a 
doubly-dispersive fading channel. This is of interest in all situations 
where spread-spectrum signals entering the channel are to be 
detected as, for example, in cell-search procedures, in handoff man- 
agement algorithms, and neighbor discovery phases. The proposed 
detection structure (referred to as MGLRT) is based on the method 
of sieves and is fully blind, meaning that it only requires knowledge 
of the spreading code of the user to be detected. We have provided 
an efficient QR-based procedure to evaluate the decision statistic 
and have also proved that the MGLRT has a bounded CFAR under 
the design assumptions. The proposed detector achieves satisfactory 
detection performance in doubly-dispersive channels, even when the 
user to be detected is only present in a fraction of the processing 
window and strong cochannel interference is present. 
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Figure 2. Pa versus SNR for a = 0.1,0.3,0.5,0.7. System parameters: 
N = 15, M = 2, L = 2, Q = 120, K = 3 and SIR= dB. 



Figure 4. Pd versus SNR for SIR= —10,0, 10 dB. System parameters: 
AT = 15, M = 2, L = 2, Q = 120, a = 0.3 and K = 3. 



